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ABSTRACT 

A nonlinear system governed by x = f(x,u) with six state variables and 
three control variables is considered in this project. A set of 
transformation from (x,u) - space to (z,v) - space is defined such that 
the linear tangent model is independent of the operating point in the 
z-space. Therefore, it is possible to design a control law satisfying 
all operating points in the transformed space. An algorithm to construct 
the above transformations and to obtain the associated linearized system 
is described in this report. 

This method is applied to a rigid body using pole placement for the control 
law. Results are verified by numerical simulation. Closed loop poles in 
x-space using traditional local linearization are compared with those pole 
placements in the z-space. 
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INTRODUCTION 


In general, there are two commonly used approaches to linearize a nonlinear 
control system: one is to linearize locally about some operating points, 

the other is to linearize globally. However, difficulties will arise when 
these approaches are applied. For the former one, it is hard to design a 
unified control law over the entire set of operating points; on the other 
hand, for the later approach, usually, the control does not permit an 
acceptable closed-loop system behavior over a wide operating range. But, if 
a nonlinear system has a tangent model independent of the operating point, 
the above stated difficulties disappear. Using this concept, Reboulet and 
Champetier [1] developed a procedure which transforms a nonlinear system 
with one control variable to another space such that the desired independent 
property is satisfied. In this project, we generalize their procedure to 
the case where a nonlinear system consists of six state variables and three 
control variables. This method is then applied to a rigid body dynamic 
model using pole placement for the control law. 


XXIII - 1 



OBJECTIVES 


The objectives of this project are: 

1. Develop a pseudo-linearization technique to analyze a nonlinear 
control system with six state variables and three control variables. 

2. Apply this technique to a rigid body dynamic model. 

3. Verify the feasibility of this technique by computer simulation. 
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BODY OF REPORT 


Concept 

Consider the class of 6-state and 3-input nonlinear systems governed by 

x = f(x,u) (1) 

where x e R 6 , u e R 3 , and f = (f^, .... fg) C : R 9 -*■ R 6 of class C 1 . The 
set of operating points is defined by 

M A {(x_,u ) : f(x ,u ) =0} 
x,u — wo o o 

and its projection in the state space is defined by 

M A{x :3u s.t. f(x ,u ) = 0}. 

X — o o o o 

In the neighborhood of an operating point (x ,u ) , the dynamic behaviour of 
the system may be considered as linear. It 2s ?hen described by 


6x = F(x ,u )Sx 4- G(x ,u )6u 
o o o o 

where 6x e R 6 , 6u e R 3 and F _A 3f/3x, G A 3f/3u. 
Our aim is to find mappings 

z ± — T^(x) , i = 1, ...» 6 


( 2 ) 


(3) 


v j = Vj (x>u) 


, i - 1, 2, 3 

(T^, . .., T^ being functionally independent and 3T, + ./3u 4 0, j = 1, 2, 3) 
such that, in the z-state space, the linear tangent Aodel is independent of 
the operating point and can be written under the controllability canonical 
form 


6 *1 

= Sz. 
4 

<5z„ 

= 6z 

2 

5 

<5z. 

= 6z, 

3 

6 


= Sv 1 

6z 5 

= 6v 2 

6 *6 

= 6v 3 


( 4 ) 
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Note that the controllability of the linear tangent model is invariant 
under the mappings (3). Hence, to obtain (4) the following condition 
must hold: 


For any (x ,u ) £ M the pair (F(x ,u ), G(x ,u )) is 
o o x,u r o o o o 

controllable. 


( 5 ) 


In practice, the solution provided by our method is effective along M , 
except at the points (x o ,u q ) for which (5) is not satisfied. 


Nuw, the variations in the z-state space 6z will be rewritten in terms of 
6x and 6u. From eqn. (1) and (3) we have 

z^ = (3T i /3x)f (x,u) , i = 1, .... 6. 

Hence, at any points (x ,u ) of M 

J v o o x,u 

6z^ = a^F(x Q ,u o )6x + ouG(x q ,u o )6u, i = 1, .... 6 


where a. are the 1-forms over M 
1 x 


a ± A d T ± 


M 


, i— 1, ..., 6. 


( 6 ) 


For convenience, reference to the operating point will be omitted. To obtain 
eqn. (4) it is easy to show the following equations must be satisfied on M 


x,u 


a.G = 0 

l 


. i - 1. 

2, 

3 

(7a) 

“i F - Vi 

» i - 1. 

2, 

3 

(7b) 

“3+1 * “i (F G) 

» i = 4, 

5, 

6 

(7c) 

with ou A d 

M 

x,u 

.3 = 7, 

8, 

9. 



With this result, the problem then becomes that of finding 1-forms a^, . . . , 
(resp. cu , (Xg, a^) satisfying eqns. (7) at any point of M (resp. M ) and 
such that there exist mappings T^(x) (i = 1, .... 6) and ^.(x,u), (f’= 7, 8, 9) 
such that 1 ^ 
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M 

x 


i = 1 


» • • * » 


6 


“i ■ d T i 


“j - d T j 


M 

x,u 


Analysis 


, j = 7, 8, 9. 


Recall that M x is a set of points described by six nonlinear equations with 
nine unknowns . ’ 

Let M be the set of operating points in the neighborhood of which the linear 

tangent model of system (1) is controllable. Since the uncontrollability of 

the tangent model (2) is equivalent to the nullity of all 5-dimension minors 

(Kalman rank condition). Therefore, the complement of ft is a closed subset 

in M which implies M is an open subset in M * 

x,u r x,u r x,u 

Furthermore, another controllability condition is given by 

rank(F-Xl G) = 6 ,\/ X . 

In particular, for X = 0 

rank(F G) = 6. 

This means that the mappings f^, i = 1, ..., 6 are functionally independent 

at the point of interest. Hence, M is a 3-dimensional submanifold of R 9 . 

Which also implies that M is a 3-dimensional submanifold of R 9 since 

dim M = dim M if M X ’£ 4> . 
x,u x,u x,u 

Let M be the projection of ft on the x-state space. In order to develop 
globai results we will suppose’ihat the following condition holds: 

TCiere exists a C 1 -diffeomorphism <|> = (<j)^, ..., <P^) t from R 6 to R 6 such 
that is given by 

<f> 4 (x) = <J> 5 (x) = <J> 6 (x) = 0 

Under this assumption, if we define 

z ± ~ » i = 1» •••» ^ 
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then the surface M x is transformed in the z-space into a 3-dimensional 
manifold D : z, = z, = z, = 0. Let a be the image of a 1-form a = d T 
under <j>, t^en 


M 

x 


a(z) A a(4> ^(z)) (3x/3z) n , z e D . 

— . U Z 

z 


(9) 


Note that if we post-multiply a = d T 



by 


(3x/3x) 


D 

z 


» 


we obtain 

D Ox/3z) , V z e D z , (10) 

z z 


a(4> (z)X3x/3z) 


D 


= d T($ (z)) 


i.e. , 


a = d 



with 


X = T ° 



(ID 


( 12 ) 


Therefore, the integrability of a over M x in the x-space can be transposed in 

terms of integrability of a over D . 

z 


Note that a only depends on z^, Z£, z^ over D z> i.e., 

6 _ 

a(z) = E a i (z^, Z 2 » z^)dz^ , z e D^. 

If we choose the 1-form a such that a^, o^, and depend on z^, Z 2 » and z^ 
only, respectively. We can define " ' " “ J 

6 

T(z) = Z / a i (z 1 , z 2 > z 3 )dz i (13) 


and we have 


6 _ 6 3 _ 

dT(z) == Z a i ( z i * z 2 » z 3 ) dz x + £ z ^[ ^ (3a i /3z .)dz . ] 
i=l i=4 j=l 3 3 

6 _ 

‘JlVV V z 3 >dz i 


a 
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Since the second term of the right hand side vanishes over D . Therefore, 
a is integrable over D and so is a over M . 

Z X 

Method 


The procedure for searching for these transformations can be decomposed into 
three steps. 

Step 1 : Find three independent directions a^, a^. perpendicular to G. 

Step 2 : From eqns. (7b) (resp. (7c)) , a^, a^» (resp. a^, ag, a^) 

can be computed. ' 

Step 3 : Integration of a^, ...,o{g (resp. a^, a g, a^) along M x 

(resp. M ) provides the desired mapping T 1 , . . . , T Q . 

X jU 1 j 


Application 

Consider a rigid body in which the angular velocities are measured along a 
body fixed axis, and with an Euler angle sequence (space - three 1-2-3) 
defining the attitudes in inertial space: 

li) + U X Id) = u 

®1 = U X + ( Vl + “z c l )s 2 /c 2 


(14) 


0 2 = U) yC;L - Vl 

0 3 = (0) y s 1 + 0) zC;L )/c 2 

where w = (u^, o) y , a^) is the angular velocity vector 

0 = (0p 0 2 , 0g) is the Euler angle vector 

u =_(u , u , u ) is the external torque acting on the rigid body 

x y z = ~ 


s^ = sin 0^ c^ = cos 0^, s 2 = sin © 2 , c 2 = cos 


and I, a 3x3 matrix, is the inertial tensor. 

Without loss of generality, principle axes are used to describe this system. 
Therefore, we can assume the matrix I is a diagonal matrix, i.e.. 
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0 


0 


I = 


0 10 

y 


ooi 




Using principle axes, system (14) become: 


03 = [u„ -(I -I )o3 03 ] — 

x 1 2 z y y z I 

x 


03 = [u -(I -I )U3 03 ] — 

y L y x z x z I 

y 


03 = [u -(I -I )03 03 ] — 

z z y x xy J I 


0, = 03 + (03 S. + 03 C-)s /c 

1 x y 1 zi-'z'z 


0_ = 03 C, - 01 s, 
2 y 1 z 1 


(15) 


0 3 = (03 yS;L + a) zC;L )/c 2 

The set of operating points is 

M x,u = {( °’ °’ °» 0 1» 6 2’ 0 3* °» °» °)It < 6 2 < 
and its projection in the state space is 

M x = {(0, 0, 0, Q v 0 2 , 8 3 )|"5 < e 2 < J) 


On M 


x,u. 


F < x o* u o ) = 


0 

0 

0 


0 

0 

0 


0 0 
0 0 
0 0 


S 1 S 2^ C 2 C 1 S 2 ^ C 2 0 0 


U 


-S, 


S l /C 2 


c l /c 2 


0 

0 


0 

0 


(16) 


(17) 
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and 


G(x o’ U o ) 


L/I 0 

x 

0 1/1 

y 

o o 



0 0 
0 0 


0 

0 

0 


J 


It is clear the system is controllable at every point (x ,u ) £ M 

o o x,u 

Now, we shall choose three independent vectors a^, cx^ all orthogonal to G. 

Let a 1 = (0, 0, 0, 1, 0, 0), a 2 = (0, 0, 0, 0, 1, 0) and a 3 = (0, 0, 0, 0, 0, 1) . 


From (7b) and (7c), we obtain 

— (1 » * C 1 S 2^ C 2* " 


a 5 = (0, c x , - S;l , 0, 0, 0) 

— (0> s^/ c 2» c j/ c 2’ • 

a 7 = (0, o, 0, 0, 0, 0, 1/I X , S;L S 2 /I y c 2 , CjSj/I^Cj), 

a 8 = (0, 0, 0, 0, 0, 0, 0, c^Iy, - s^IJ, 

a Q = (0, 0, 0, 0, 0, 0, 0, sjl c„, cjl c_). 

y l y 2 i z 2 

After integrating a., i=l, .... 6, (resp. a , a_, a.) over M (resp. M ), 
the desired transformations are: X x,u 
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W x + ( Vl + Vl )s 2 /c 2 

<Vl " Vl 3 = ®2 

• 

(w s. + a ) c.)/c 0 = 0, 
y 1 z 1 2 3 


v, + v,, 
la lb 


v 2a + v 2b 


v 3a + v 3b 


= u z n x + (u y s l 7l y + u z C l^ I z )s 2 /c 2* 

- {[(I Z - VVA + (I x - I y>Vl /l zK S 2 /c 2 + 

(I y " I z )u) y w z /I x } + (u) y C l ' “z S l ) K S 2 /c 2 + (w y S l + X 

(1 + 2s 2 2 /c 2 2 ) , 

= u c./l - u s,/l , 
y 1 y z 1 z 

= [(I - I )w c,/I + (I - I )u> s n /l ]o) - 
z x 7 z 1 y x y y 1 z x 

(o) y S 1 + W z C 1 )[( i ) x + (UySj^ + W z C 1> S 2 / C 2 3 » 

= CUySi/i y + u^/y/cj, 

= [(I - I >03 S./I + (I - I )u) C./I ]U) /c + 

z x 7 z 1 y x y 7 y 1 zz J x z 

(w c-j^ - W z s 1 )[(i3 x + (to s 1 + 0)^X1 + l/c 2 )s 2 /c 2 ]/c 2 
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Note that each transformed input v. (i^ =1, 2, 3) is partitioned into 
v^ a and v., terms. The v^ a term arises directly from integration of the 
corresponding ct^ + . vector. The v., term can be identified with the 
nonlinear terms wftich appear in the expression of According to [1], 

these high order terms may be neglected in the neighborhood of the operating 
point, in which case the equations of motion are said to be psuedo-linearized. 
However, in this particular case, since the operating points require w = 0, 
the transformation from v. can be augmented to include the v., term without 

f J j_ 1 — i. I _ • * i £ rrru ^ 


affecting the requirement that dv^ 


= a., i = i+6. Therefore, as shown 
M j J 

x,u 

before, the equation of motion for a rigid body become extirely linear when 
expressed in terms of Euler angles and Euler angular rates. 

Taking account of (18), (15) can be rewritten as a linear systems: 


Z 1 ~ V 




z„ = z 


6’ 


z 4 = v l* 


z 5 V 2 * 


z 6 “ v 3* 


or, in a matrix form: 


V. Z 6 J 


0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


1 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 


N 

0 


<■ N 
z 






0 

0 

0 

0 

1 

0 


V 1 


(19) 


) 


Furthermore, the inverse transformations also can be found as: 
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- Z 5 S 1 + Z 6 C 1 C 2 


(z 5 + U) z s 1 )/c 2 


2 4 - <M y s l + “z c l> s 2 /c 2 


V <v l ■ ’lb) - (v 3 - v 3b )s 2 ] 


I y t(v 2 - v 2b )c l + (v 3 " v 3b )c 2 s l 1 


V c l c 2 (v 3 ' v 3b> - s l (v 2 - v 2b) ] 
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The above transformations were verified by numerical integration in both 
the x-space and z-space for a prescribed input u. Trajectories obtained 
by integrating in x-space were virtually identical to those obtained by 
integrating in z-space then transforming to x-space via equation ( 20 ). 
Figure 1 shows the effect of psuedo-linearization at an operating point 
of ©i = ©2 = 63 = ir / 4 with a unit step input for u. The solid lines 
are trajectories obtained by directly integrating equation (15), while 
the dashed curves represent the transformation from z-space to x-space 
neglecting the v terms. 


For control studies, a constant gain, full state feedback control v = -Kz 
was employed. Working in the z-space with z = Az + Bv, the following 
gain matrix K allows for arbitrary, complex conjugate pole placement, 


K = 


0 )^ 0 


0 a). 


0 0 (jo- 
's. 3 


o 2 e l w L 


o 

0 


2 ? 2“2 


0 1 


( 21 ) 


2 S“3 


where and prescribe the closed loop frequency and damping ratio of 
the ith axis. The actual control to be used in x-space comes from equations 
(18), ( 20 ) and ( 21 ), where both u) and© are assumed to be available as plant 
measurements. Figures 2, 3 show the system response to the arbitrary initial 
conditions = ©2 = O 3 = tt/ 4 and w x = d)y = to z =0.5 using control frequencies 
and dampings of U)^ = = 1 .0 and ~ *>2 = ^3 = 0*707. Once again, 

the solid lines represent tne response using the exact transformation of 
equation ( 20 ) while the dashed lines show the effect of psuedo-linearization. 
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OPEN LOOP EFFECTS OF PSUEDO-IINEARIZATION 



TIMECSEC) 
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OPEN LOOP EFFECTS OF PSUEDO-LINEARIZATION 



TINECSEC) 
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u 3> m zc - j* ru3>Hmn'* Mi>Hrnn 


CLOSED LOOP EFFECTS OF PSUEDO-LINEARIZATION 



TIHECSEC ) 
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UN 


CLOSED LOOP CONTROL IN THE 2 - PLANE 



TINE(SEC) 
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CLOSED LOOP CONTROL IN THE 2 - PLANE 



TIME (SEC) 
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CONCLUSION AND RECOMMENDATION 


A global linearization method for a nonlinear control system with 6 states 
and 3 inputs is successfully developed in this project. This method is an 
extension of a pseudo-linearization technique proposed by Reboulet and 
Champetier [1] in which only one input is considered. After applying this 
method to a rigid body in which the angular velocities are measured along 
a body fixed axis, and with an Euler angle sequence defining the altitudes 
in inertial space, the original non-linear system is described by a linear 
system in a transformed space. Therefore, a global control law can be 
established accordingly. 

For further application of this method, we will turn our attention to the 
case of control systems of multiple bodies. Applying this method, Mr. Sharkey 
has obtained very positive numerical results on a two-body model using pole 
placement for a control law. We suggest a mathematical research to support 
this new application which is necessary and should start immediately after 
this program to pave the road for the next summer’s project. 
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